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Abstract 

In contrast to other approaches to (2+l)-dimensional quantum gravity, the 
Wheeler-DeWitt equation appears to be too complicated to solve explicitly, 
even for simple spacetime topologies. Nevertheless, it is possible to obtain 
a good deal of information about solutions and their interpretation. In 
particular, strong evidence is presented that Wheeler-DeWitt quantization 
is not equivalent to reduced phase space quantization. 
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In the continuing search for a realistic quantum theory of gravity, it has often proven use- 
ful to explore simpler models that share the basic conceptual features of general relativity. 
Gravity in 2+1 dimensions is one such model, providing a fully diffeomorphism-invariant 
theory of spacetime geometry that nevertheless avoids many of the technical difficulties of 
realistic (3+l)-dimensional gravity. 

The goal of this paper is to explore the ramifications of one popular approach to quan- 
tum gravity, the Wheeler-DeWitt equation, in this simple setting. We shall see that even 
with the simplifications of 2+1 dimensions, Wheeler-DeWitt quantization is considerably 
more complicated than one might guess, and that at least in its simplest interpretations, 
it is not equivalent to other known approaches to quantization. 



1. Canonical Gravity in 2+1 Dimensions 

Before trying to construct a quantum theory, let us briefly review the canonical formal- 
ism for (2+l)-dimensional gravity, as described by Moncrief JTJ and Hosoya and Nakao JJ. 
We begin with a spacetime with the topology RxE, where £ is a closed surface. The 
standard (2+l)-dimensional ADM variables are then a spatial metric g^ and its canonical 
momentum 7r y . By a classical result of Riemann surface theory, any metric on £ can be 
written in the form 

9ii = e 2X f*9ij, (1-1) 
where / is a spatial diffeomorphism and is a "standard" metric of constant curvature 1 
(if £ is a sphere), (if S is a torus), or —1 (if S is a surface of genus g > 1). For a surface 
of genus g > 1, the standard metrics g^ comprise a (6g — 6)-dimensional family; for the 
torus, the family is two-dimensional, and we can choose 

ds 2 = T 2 - 1 \dx + Tdy\ 2 , (1.2) 

where x and y are angular coordinates with period 1 and r = Ti+ir 2 is a complex parameter, 
the modulus. Corresponding to the decomposition ( p. . 1|) , the momentum tt^ can be written 

as 

7r« = e" 2A v^ (p ij + i^/vf + VY J + V j Y* - g ij V k Y h \ , (1.3) 

where V, is the covariant derivative for the connection compatible with g^, indices are now 
raised and lowered with g^, and p 1 ^ is a transverse traceless tensor with respect to V, (in 
the language of Riemann surface theory, a holomorphic quadratic differential). Roughly 
speaking, is canonically conjugate to g^, tx to A, and Y l to /. More precisely, if we 
consider a cotangent vector 5gij in the space of metrics, 

6 gij = VMj + + 2e 2X g lJ 5X + e 2A <%, (1.4) 

the symplectic structure can be read off from the expression 



S gij 8ir ij = |VP gij Sp ij + 25 X 5tt-2 5? 



e- 2X ^j (A + ~r) 5Y t + I Vi (V 2A 5tt 
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(Once again — and for the remainder of this paper — indices are raised and lowered with 
gij.) The last term in ( |1.5|) tells us that £ l and Y{ are not quite a conjugate pair, a fact 
that will cause us considerable grief in the next section. 

In terms of the decompositions fll.lQ an d ( |1.4|) , the constraints of canonical gravity 
become relatively simple. The momentum constraints are [|J 

V?(A + ±R)Y t + \e 2X Vi (e^vr) = (1.6) 
- note the similarity to the last term of (|1.5| ) — while the Hamiltonian constraint becomes 



H = ~\^- 2X K 2 + ^e- 2A [p ij Pij + 2 Pij {PYf + [PY) l3 {PYf] + 2^ 
where 



AX--R 
2 



= 0, 

;i-7) 



{PY) {j = ViY 3 + VjYi - g i3 V k Y k . (1.8) 

Moncrief has shown that in the classical theory, these constraints fix A uniquely as a function 
of cjij and p % \ thus determining a finite-dimensional reduced phase for (2+l)-dimensional 
gravity. A description of dynamics on this reduced phase space depends on a choice of 
time slicing; for a foliation by surfaces of constant mean curvature Tri^ = T, it may be 
shown fllLBI that the Hamiltonian is 



H = Jd 2 xe 2X , (1.9) 
which can in turn be written as a (complicated) function of and . 



2. Quantization and the Wheeler-DeWitt Equation 

To quantize this system, we may either solve the constraints classically and quantize 
the resulting reduced phase space, or else quantize first and impose the constraints as 
conditions to determine physical states. The first alternative is discussed in reference ||. 
If S has the topology of a torus, the constraints can be solved explicitly; in particular, in 
the York time slicing, 7re~ 2A / y/g = T, the momentum constraints imply that Yj vanishes, 
and the Hamiltonian constraint has as its solution 

e 2A = T- 1 (2 Ft ,^) 1/2 . (2.1) 

The system is described by the effective Hamiltonian ( |1.9[ ), and it is an easy exercise to 
find the corresponding Schrodinger equation: 

«'H = T- l A ^, (2.2) 
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where 



A = -r 2 2 



drr 



+ 



a 2 ' 
<9r 2 2 



(2.3) 



is the Laplacian on the torus moduli space with the standard Poincare metric dr<if/r 2 2 . 
The square root in (|2.2j ) is ambiguous, of course; it is customary to take it to be the unique 
positive square root, corresponding physically to an expanding universe. 

In Wheeler-DeWitt quantization, on the other hand, the Hamiltonian constraint is 
not solved classically, but is instead imposed as a Klein-Gordon-like equation of motion 
constraining physical states. In this approach, we no longer have the luxury of choosing 
a time slicing, and in particular cannot require that TrK be constant on each slice. As a 
consequence, the momentum constraints no longer require Yj to vanish, but imply instead 
that 

1 ° x v,(, 



Yt = --A- 1 
2 



V9 



e 7T 



(2.4) 



To obtain the Wheeler-DeWitt equation, we are instructed to insert this expression into 
the Hamiltonian constraint and to make the further substitutions 



7T = — 



% 5 
26X'' 



V 
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(2.5) 



For the torus, we find that 



l^ 2X ^ + l^ 2X ^o + 2AA - 2e- 2 ^[n]AY t [n] 



(2.6) 



+ 2e" 2A V 



(2p ij + V^tt] + V j Y l [n] - g lj V k Y k [7i}) Yj[v] 



'*[A,r] = 0, 



with Yi[ir] given by Q2.4p . The operator ordering in the first term of (|2.6| ) is the simplest 
for which the Hamiltonian constraint is at least formally Hermitian with respect to the 
inner product defined by the "measure" [dX\. This is a natural choice, but it is by no 
means unique; we shall return briefly to this issue in section 5. As promised, the Yj 
dependence makes equation ( |2.6| ) complicated and nonlocal, reflecting the fact that the 
symplectic structure mixes Y^ and 7r — only in the York time slicing do the momentum 
and Hamiltonian constraints disentangle. 



3. Naive Schrodinger Interpretation 

It is evident from equation (|2.6| ) that explicit solutions of the Wheeler-DeWitt equation 
will be difficult to find, even in this simple model. The fundamental simplifying feature 
of 2+1 dimensions is that the physical phase space, parametrized by and p l \ is finite 
dimensional, effectively reducing quantum field theory to quantum mechanics. This sim- 
plification is evident in the Pijp 1 ^ term of the Hamiltonian constraint, which reduces to a 
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Laplacian on the two-dimensional moduli space of the torus. But unlike other approaches 
to quantization, the Wheeler-DeWitt equation retains a complicated dependence on the 
scale factor A, reflecting the fact that the time slicing has been left arbitrary. 

Even without an exact solution to the Wheeler-DeWitt equation, we can try to extract 
some useful information. To do so, however, we must first decide on an interpretation of 
the "wave functions" \I/[A, r]. This is a highly nontrivial problem: solutions of ( |2.6| ) do not 
come with a ready-made inner product, and the choice of something like a Hilbert space 
structure is necessary to give wave functions a quantum mechanical interpretation. 

One natural guess, the "naive Schrodinger interpretation" @], is that |^[A, r]| 2 is simply 
the relative probability of finding a spatial slice of the universe with scale factor A and 
modulus r. To compare Wheeler-DeWitt and reduced phase space quantization, it is 
more useful to consider the relative probability of finding a slice with extrinsic curvature 
K = e~ 2X 7T and modulus r; this can be obtained from a functional Fourier transform 

%[K, t} = J [dX] e ^ d2xKe2X m[\ t\. (3.1) 

In particular, the wave function on a York time slice K = T is 

ijj(T, t)= J [dX] e iT I d2xe2X ^[X, r], (3.2) 

giving the relative probability of finding a geometry with modulus r on a slice of constant 
mean curvature T. 

We can now ask whether this wave function is equivalent to the wave function (|2.2| ) 
obtained by fixing the same time slicing before quantization. Observe first that by ( |2.4| ). 

J [dX] eVf^YifrWX, r] = 0, (3.3) 

since functional integration by parts can be used to move the action of it to the exponential. 
The functional Fourier transform of the Wheeler-DeWitt equation fl2.6|) then reduces to 

J [dX] e iT I d2xe2X J d 2 x N(x) {-T 2 e 2A + e~ 2A A + 4AA} #[A, r] = 0, (3.4) 

where N(x) is an arbitrary function. The expression in brackets is essentially the classi- 
cal Hamiltonian constraint for the torus, with PijP^ replaced by A . Now, however, the 
constraint need not be obeyed everywhere, but only inside a functional integral. In partic- 
ular, there is no reason to suppose that the wave function \l/[A,r] has its support only on 
solutions of this constraint. 

To obtain further information about ip(T, r), we can set N = 1 in ( |3.4| ). We find that 

^T 2 ^(T,r) = -A /[dA] [e^I^ Jd 2 x e~ 2A ] ¥[A, r], (3.5) 
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or, differentiating again, 



gfT — + A 



ip(T, r) = J [dX] | e lT /^ e2A [l - J d 2 x e 2X J d 2 x e' 2X ] } A *[A, r 



(3.6) 

Up to minor operator ordering ambiguities, the left-hand side of this expression is the square 
of the Schrodinger equation ( |2.2| ) obtained from solving the constraints before quantizing, 
and is equivalent to the "gauge-fixed Wheeler-DeWitt equation" obtained by Hosoya and 
Nakao []2| and Martinec ||. The right-hand side thus measures the departure of the full 
Wheeler-DeWitt equation from these alternative approaches to quantization. 

In particular, the right-hand side of ( |3.6p vanishes when A is spatially constant, as it 
is for any classical solution of the constraints. As we have seen above, however, there is 
no reason to expect the wave function \I/[A, r] to have its support only on such solutions. 
Moreover, the right-hand sides of ( |3.5| ) and ( |3.6| ) depend not only on ip(T, r), but on the 
full Wheeler-DeWitt wave function \I/[A, r], so our would-be Schrodinger equation must 
be interpreted as one of an infinite family of equations — somewhat analogous to the 
Schwinger-Dyson equations in quantum field theory — for functions 

MT, t)= f [d\] U T I^ X f d 2 xe 2nX ) *[A, r]. (3.7) 



While the particular form of equation (|3.6|) depends on a choice of operator ordering in the 
Hamiltonian constraint, I have been unable to find any alternative ordering in which the 
Wheeler-DeWitt equation reduces to a single Schrodinger- or Klein-Gordon-like equation 
for ip(T, r). 

The details of this argument depend on a specific model, but the general features seem 
likely to extend to realistic (3+l)-dimensional gravity. In particular, the trace of the 
extrinsic curvature in 3+1 dimensions is canonically conjugate to y/g, so the wave function 
at fixed K should take the form 

^ ( 3+i)[T,&] ~ J [d(^)}e iTV ^ (3+1) [^,~grj}, (3.8) 

where V = J d 3 Xy/g is the volume of a spatial slice and <?jj symbolizes the remaining metric 
components. If we take two time derivatives of ^^3+1), the integrand will be multiplied by 
a nonlocal factor ^ 2 [y / ^]. But just as in our (2+l)-dimensional model, the action of the 
local part of the Hamiltonian constraint can generate at most a single spatial integral in 
the integrand. Unless the nonlocal terms in the Hamiltonian constraint conspire to remove 
this mismatch in the number of integrations — as they fail to do in (2+1) dimensions — it 
would seem difficult to obtain a simple Klein-Gordon-like equation for the wave function at 
constant K. The inequivalence of Dirac and reduced phase space quantization in a theory 
with quadratic constraints has been widely studied — see, for example, ~ but I 

believe this particular feature is new. 
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4. Gauge-Fixing the Inner Product 



The absence of a simple Schrodinger or Klein-Gordon equation for ip(T, r) should not 
in itself be viewed as a reason for rejecting Wheeler-DeWitt quantization. There are many 
possible quantum theories of gravity in 2+1 dimensions [P0||, and we do not yet know how 



to choose among them. The naive Schrodinger interpretation has a somewhat more serious 
problem, however — there is no reason to expect the inner product (ip{T, r)|0(T, r)) of two 
states to be conserved, so the relationship to ordinary quantum mechanics is problematic. 
This difficulty can be resolved, at least formally, by appealing to a suggestion by Woodard 
|Tl"|j for redefining the inner product. 

Woodard's proposal is that all invariances of the action, including those generated by 
the Hamiltonian constraint T~C, should be gauge-fixed in the inner product of two Wheeler- 
DeWitt wave functions. For transformations generated by the diffeomorphism constraints, 
this procedure is uncontroversial: such transformations take a metric on a spatial slice £ 
to a physically equivalent metric, and must clearly be factored out to avoid overcounting. 
The proper treatment of the group of transformations generated by the Hamiltonian con- 
straint, on the other hand, is less obvious. Such transformations mix metrics on X with the 
corresponding canonical momenta — in the language of geometric quantization, they do 
not preserve the polarization — and they have no simple interpretation as transformations 
of the argument of the wave function. 

One argument in favor of Woodard's approach comes from examining wave functions 
determined by path integration. Consider two Hartle-Hawking wave functions ^i[g] and 
^2[g] on a surface X, obtained from path integrals on manifolds Mi and Mi with boundaries 
dM\ « DM2 ~ X. On Mi, say, the action is invariant under the transformations 

Sgij = [eH,gij} 

57c ij = [eH,7r lj ] (4.1) 

only when e = on X ||12|| , so only this subset of transformations will be gauge-fixed in 
the path integral for \I/i[g]. A similar statement holds for ^[f?]- The inner product, on the 
other hand, can be naturally defined as the path integral over the new manifold M obtained 
by attaching Mi and M 2 along X. In this path integral, however, the transformations ( |4.1| ) 
must be gauge- fixed no matter what the value of e on S; this additional gauge- fixing leads 
directly to Woodard's formalism. 

To apply this approach to our (2+l)-dimensional model, we begin with the inner product 



= J[d\] J -L**[\,t]®[\,t}, (4.2) 
gauge-fixed to the York time slicing 

X = n- Te 2X = 0. (4.3) 
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The path integral over A can be evaluated a la Faddeev and Popov by inserting a factor of 
unity in the form 

1= [[de]v*[\*,T e ]5[x e M* € ,T% (4.4) 



where A e , for instance, denotes the transformed value of A, 

A £ = A - i [eH, A] . (4.5) 

A simple computation shows that 



z/[A, r] = det 1/2 



1 „o o \ 1 _ x 





2 



(4.6) 



up to operator ordering ambiguities coming from the dependence of x on the canonical 
momentum 7r. Note that unlike the more familiar Faddeev-Popov determinants of gauge 
theory, z/[A,r] is an operator, depending explicitly on the Laplacian A on moduli space. 

Finding the gauge-fixed form of the inner product ( |4.2| ) is now straightforward: it is 
easy to check that 

<*|*> = J ^r(T,T)t(T,T) (4.7) 

with 

■0(T, r) = / [dX)u[X, r] e iT f d2xe2X ^{\, r). (4.8) 



The wave function tp(T, r) is almost identical to the naive Schrodinger wave function ^(T, r) 
of the previous section, differing only by the presence of the determinant z/[A,r]. Now, 
however, the inner product (|4.7| ) is guaranteed to be conserved: the parameter T merely 
labels an arbitrary choice in the gauge- fixing of the original path integral ( |4.2j ), and the 
final result must be independent of any such choice. 

It is natural to ask how the the added factor z/[A, r] affects the arguments of the previous 



section. Unfortunately, this is a hard question: equation (|4.6j ) is complicated enough that 
I have been unable to find a useful explicit expression for the determinant. Observe, for 
instance, that the nonlocal terms in the Wheeler-DeWitt equation ( |2.6| ) can no longer be 
neglected. Moreover, u[X, r] depends on the modulus r, and does not commute with the 
Laplacian A in the Hamiltonian constraint; to obtain the analog of the would-be Klein- 
Gordon equation (|3.6|) , we would have to know this dependence explicitly. 

Nevertheless, we can at least obtain some qualitative information about the determi- 
nant ^[A, r] without too much difficulty. Recall that the naive Schrodinger wave function 
ip(T, t) failed to obey a simple Klein-Gordon equation because the Wheeler-DeWitt func- 
tional \I/[A,t] had support away from solutions of the classical constraints. It is therefore 
interesting to note that ^[A,r] is itself peaked around solutions of the constraints. Indeed, 
let us write 

D = —A + Ve 2A + ^ 2A A = D + (V - V ) (4.9) 
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so 

V 



with 

V = I T 2 e 2A + I e - 2A A , V = Jd 2 x V , D = -A + V . (4.10) 
Note that Vo is always positive, so Do is invertible. It is not hard to show that 

{oo (_-i\n+l ~\ 

But D has eigenfunctions and eigenvalues 

Ait 2 

\ mn ) = e^irnx+ny) ? ^ = | n _ mr |2 + ^ 

[A, r] = det^pol exp £ -^i|(mn |V - V | W )| 2 + •••}, (4.13) 

^ ^ 'inn ^pq ) 

which is clearly peaked at (mn\V — Vo\pq) = 0, i.e., at spatially constant values of A. 
Moreover, for A constant, the derivative 

Irdet^l^l = (r 2 e 2X - e- 2A A ) • J- det 1 / 2 ^ (4.14) 
cm v ' oVo 

vanishes precisely at the solutions ( [2.1|) of the Hamiltonian constraint, now interpreted as 
an operator equation on moduli space. One can further check that in any regularization in 
which TtDq 1 is positive, this extremum is in fact a maximum. 

Woodard's proposal thus leads to wave functions at fixed TrK that receive their main 
contribution from configurations that nearly satisfy the classical constraints. In particular, 
we see from Q4.13 ) that deviations from a constant scale factor are exponentially suppressed. 

1/2 

Note that for A near its classical value, V ~ TA , so the eigenvalues £ mn will be smallest 
- and the suppression strongest — at small values of T, i.e., late times. [] The analog of the 
would-be Klein-Gordon equation ( |3.6|) for Woodard's wave functions should therefore be 
closer to the true Klein-Gordon equation of reduced phase space quantization, especially 
at late times. 



5. Next Steps 

Can the results of the last section be made more precise? Perhaps. The two key prob- 
lems are to understand the determinant u[X, r] and the operator ordering in the Hamiltonian 
constraint Ti,. For any fixed A, u[X, r] is a modular invariant function on the torus moduli 

*In an expanding (2+l)-dimensional universe, the value of TrK on a constant mean curvature slice 
decreases from infinity at the initial singularity to zero in the far future. 
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space, and as such can be expanded in terms of eigenfunctions of the Laplacian Ao- For 
instance, it is not hard to show that the zeta function for the operator D is 



U M = v - + (40- 1 f^T) (-M E ^ * + *>■ P' 1 ' 



where E(s, k) is the Eisenstein series [p~3] , pT4 



E(r,a)= E / , 2 r (5-2) 

771 + TT1 



(The prime means that the point m = n = is omitted from the sum.) Such series have been 
studied quite extensively, and equation (|5.2|) may yield some useful information. Elizalde 
has pointed out that the sum may be further simplified by an analog of the Chowla- 
Selberg relation [[Kfl. Similarly, the exponent in ( [4.1 1| ) can be expressed in terms of another 
Eisenstein series fl~4}| , 



E*, v (t,s)= e 2 ^™ +ny \ T% .. . (5.3) 

' — \m + rn\ 



These quantities are complicated, but perhaps not completely intractable. 

As for operator ordering, one possible approach is to search for an ordering such that 
[H(x),7i(x > )} = 0. This is a highly nontrivial task, however, even without the nonlocal 
term in ( |2.6| ). For example, the commutator of Ao and A A is a complicated expression 
involving single derivatives with respect to the modulus r, and there seems to be no simple 
ordering which eliminates these terms. Once again, the Wheeler-DeWitt equation proves 
to be more complex than one might have anticipated from other approaches to (2+1)- 
dimensional gravity. 

Using this model to investigate other interpretations of the Wheeler-DeWitt equation 



would also be of considerable interest. Wald [16j has recently suggested an intriguing ap 



proach to defining a Hilbert space structure, based on the analogy between the Hamiltonian 
formulation of general relativity and that of the relativistic particle. This approach requires 
a solution of the Cauchy problem for the Wheeler-DeWitt equation (|2.6|), however, a task 
made quite difficult by the presence of nonlocal terms in the Hamiltonian constraint. 
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